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Electrical and thermal characteristics of lithium-ion battery packs in electric vehicles in different operat¬ 
ing conditions are important in order to design the battery pack thermal management system. In this 
work, electrical and thermal behaviors of different size of LiFeP0 4 cylindrical cells are investigated under 
various operating conditions. The simulation results show good agreement with the experimental data 
under various operating modes. Due to the large thermal resistance of layered active material in a Li- 
ion cell, the temperature difference in the radial direction is significantly correlated with a diameter of 
cell and I t - rates. Compared with natural convection, strong forced convection will reduce the temperature 
uniformity in the cell and accelerate the thermal aging rate. Lastly, integration issues of the cells into a 
battery pack are discussed from mechanical, electrical, thermal, control and monitoring, manufacturing 
and maintenance aspects. These issues could impact the performance, cost, driving range and life cycle 
of the battery pack in electric vehicles. 

© 2014 Elsevier Ltd. Ah rights reserved. 


1. Introduction 

Li-ion batteries are rechargeable batteries which are gaining 
popularity for Electric Vehicles (EVs) and Hybrid Electric Vehicles 
(HEVs) applications. Li-ion batteries with high energy density, 
low maintenance, less toxic, good cycle life and capable of accept¬ 
ing high charging rate is suitable for energy storage system in EVs 
and HEVs. Moreover, Li-ion batteries have no memory effect, do 
not required periodic deliberate full discharge and the self¬ 
discharge rate is less than half of Nickel Cadmium (NiCd) and 
Nickel Metal Hydride (NiMH) batteries. Hence, Li-ion batteries 
have been extensively investigated for potential applications in 
EVs and HEVs to replace NiMH and lead acid battery which have 
low energy density and depth of discharge. Recently, a series of fire 
accidents involving Li-ion battery pack in EVs indicated that there 
are still many challenges to be overcome, especially thermal issues 
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of Li-ion battery. Li-ion battery needs to operate between 
25 °C-40 °C to maximize its performance, cycle life and reduce 
capacity fading caused by thermal aging [1 ]. Cycle life and capacity 
of the Li-ion battery are inverse proportional to the temperature of 
the cell [2,3]. In order to satisfy the power requirements for specific 
devices, Li-ion batteries are electrically connected in series and or 
parallel to form a pack. Therefore, the uniformity of the cell tem¬ 
perature in the large pack should be maintained as homogeneous 
as possible (3-5 °C) to ensure comparable power performance, 
effective cell balancing and charge acceptance during regenerative 
braking [4]. Therefore, an accurate battery model is needed in 
order to predict the temperature of the battery pack when the 
physical battery testing facility is not available. Battery model 
has several advantages such as estimation of electrical and thermal 
behavior under different driving conditions in a short time, reduce 
the product life cycle cost of the battery pack and testing equip¬ 
ment cost [5]. The battery model is particularly useful for battery 
pack thermal management system and battery pack control and 
monitoring design and planning. 
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Nomenclature 



A 

exponential voltage, V 

R c 

terminal contact resistance, Cl 

A s 

external surface area of the cell, m 2 

T 

temperature of battery, K 

A m f 

cross section flow area for cooling air per module, m 2 

To o 

free stream temperature of air, K 

A ms 

total module surface area exposed to cooling air, m 2 

t 

time, s 

B 

exponential zone time constant, (Ah) -1 

Vbatt 

battery voltage, V 

CFM 

Cubic feet per minute 

P 

density of the battery, kg m -3 

Cp 

specific heat capacity of the battery, J kg -1 K -1 

Pair 

density of the air, kg m -3 

Qjp 

theoretical capacity of the battery, Ah 

Gsb 

Stefan-Boltzmann constant, W m -2 I< -4 

E 

emissivity 



Eo 

battery constant voltage, V 

Subscripts 

h 

convective heat transfer coefficient, W m -2 K -1 

exp 

exponential 

h 

discharge current in amperes during 1 h discharge 

nom 

nominal 

i 

battery current, A 

full 

fully charged 

i 

filtered current, A 

max 

maximum 

I< 

polarization constant, V(Ah) -1 or polarization 
resistance, Cl 

forced 

forced convection 

2 

m 

conductive heat transfer coefficient 
mass flow rate of air, kgs -1 

Superscript 

k Peukert constant 

Q 

battery capacity, Ah 



R-int 

internal resistance, Cl 




Several Li-ion battery models have been proposed to predict the 
charging and discharging behavior, state of charge (SOC), chemical 
reaction, Li-ion distribution, current density distribution, heat 
generation, temperature distribution, etc. Li-ion battery models 
can be distinguished into two major groups, mechanistic models 
- electrochemical, thermal models and empirical models- equiva¬ 
lent circuit models. Electrochemical modeling [6-10] used a cou¬ 
pled time variant spatial partial differential equations to model 
the electrochemical reaction of the battery. Although electrochem¬ 
ical models can accurately predict the aging and thermal behavior 
of Li-ion battery, the equations itself are complex and require 
extensive computational resources [11]. Moreover, most of the 
studies only presented results of numerical simulations and did 
not validate with experimental work 12-15]. 

The electrical model [16-27] or equivalent electric circuit 
consists of the voltage source, resistors and capacitors or set of 
empirical equations to represent the electrical behavior of a bat¬ 
tery. The open circuit voltage of the battery is defined as a function 
of the SOC. The parameters required for the model are extracted 
from the Hybrid Pulse Power Characterization test (HPPC), pulse 
discharging or low J r rate of constant current discharging curves. 
These models are easy to use and required shorter computational 
times and have an advantage in large system simulations, espe¬ 
cially in analyzing dynamic and thermal behavior of the battery 
under different driving cycles. 

Lithium iron phosphate (LFP) battery have high thermal 
runaway temperature (270 °C), non-toxic, do not release oxygen 
at elevated temperatures and high cycle life (1000-2000) as com¬ 
pared to LiCo0 2 , LiMn0 4 and LiNiMnCo0 2 has make it become an 
attractive solution for future EVs battery pack [28]. Although there 
are some studies have been carried out to investigate the thermal 
behavior of Li-ion battery under different driving cycles, modeling 
work on LFP cell are rare [20,29-31]. 

Li-ion batteries can be constructed into various shapes and sizes 
for different applications: spiral wound to form cylindrical cells or 
stacked plates to form pouch cells. The fabrication technology of 
spiral wound design of the Li-ion battery is more mature and 
commonly found in the market. Besides, spiral wound design has 
several advantages as compared to prismatic and pouch cell design 
such as: easy to fabricate, high energy density, mechanical stabil¬ 
ity, incorporated safety vent and not prone to swelling during 


operation [32]. However, small specific area to volume of a cylin¬ 
drical battery could lead to the development of a large temperature 
gradient in the cell and thermal aging issues, heat is retained in the 
cell and hot spots are formed in the center of the cell. Large tem¬ 
perature gradient developed across the cell at large load could 
cause capacity fading [33]. The temperature gradient in the cell 
is significantly correlated with the diameter of the cell and J r rates. 
This implies that thermal analysis of different sizes of cylindrical 
cell is necessary to provide a guideline for determining a suitable 
size to integrate into EVs or HEVs battery pack. 

In the present work, an electrical-thermal battery model was 
developed to investigate the electrical performance and thermal 
behavior of two different sizes of spiral wound LFP cell. The 
simulation results are confirmed by experimental data. The devel¬ 
opment of internal temperature under different cooling conditions 
was predicted by the battery model. Next, the validated single cell 
model is then extended to the whole battery pack to examine the 
thermal response of the cell under US06 Supplemental Federal Test 
Procedure (SFTP) driving cycle. Furthermore, the integration issues 
of the cell into the battery pack are discussed from different points 
of view, such as mechanical, electrical, thermal, manufacturing and 
maintenance. This study will serve as a basic guideline for cell 
thermal management system design and integration of cell for 
EVs and HEVs battery pack. 

2. Model development 

2 A. Battery model 

A battery model is needed to describe the correlation between 
the input parameters to the Li-ion battery such as current or power 
and outputs of the battery model like voltage, SOC for constant cur¬ 
rent charging/discharging and dynamic conditions. The modified 
Shepherd model was used to model the voltage dynamics of the 
LFP cell in this study [34-36]. Eqs. (1) and (2) are used to model 
the charging and discharging characteristics of the battery respec¬ 
tively [35,36]. In this model, the internal resistance of the cell was 
assumed to be constant and does not change with SOC and temper¬ 
ature. Moreover, the battery model was treated as independent of 
the temperature and the model parameters for charging and dis¬ 
charging are identical and the hysteresis phenomenon was not 
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Fig. 1. Typical discharge characteristic of Li-ion battery [35,36]. 


modeled in this study. The battery model parameters E 0 , I<, A and B 
were extracted from the battery discharge curve of 0.2 I t as in Fig. 1 
in according to the procedures proposed by Tremblay et al. [35,36]. 
Filter current is used to model the slow dynamic behavior of cell 
voltage and solve the algebraic loop problem in Simulink [35]. 

Charging (f < 0) 


+ Aexp(-B • it) (1) 

Discharging (f > 0) 

V bat t =E 0 - R int i-K 3L. ■ f - K (^) it + Aexp(-B ■ it) (2) 


The voltage of the cell at fully charged state is defined in Eq. (3) 

[35,36]. 

Vfull = E 0 - Rint i +A (3) 

The voltage at the exponential section is defined in Eq. (4) 

[35,36]. 

V exp = Eo ~ I( Q_Q (Qexp + 0 — Rint ' * + Aexp Qexp^J (4) 

The voltage of the cell at the nominal zone is defined in Eq. (5) 

[35,36]. 

Vnom — Eq ~ K — — • (Q.nom + l) — Rint ' 1 

Lc — ^Lnom 


+ Aexp ( • Q.nom 


( 5 ) 


2.2. Thermal model 

In this study, the physical properties of the battery are assumed 
uniform and not affected by temperature. General energy balance 
as in Eq. (6) proposed by Bernardi et al. was used to identify the 
total heat generated in the cell [37]. The source of heat generated 
in the cell comprises joule heat and reversible heat. The reversible 
heat term is calculated using the relation of dU/dT with the SOC 
proposed by Forgez et al. [38]. In addition, contact resistance has 
been added in the equation to improve the prediction of the heat 
generation. 


The general energy balance equation for the battery thermal 
model is given by: 

r dT , n ,i , 

PCp dt — ^ g r 2 + 0) 

With boundary conditions at the outer surface of the cell (r = R) 
is defined by 

= h(T-T 00 )+Ea sb (T 4 -T*J (8) 

r=R 



2.3. Battery pack thermal model 

Modeling of EV is based on Hyundai Trajet and the battery pack 
is designed using 38,120 cells to provide similar power capacity as 
in previous studies [28]. Axial air cooling architecture was used to 
investigate the thermal response of the EV battery pack [28]. The 
battery pack comprised 28 modules and each module is con¬ 
structed using 28 pieces of 38,120 cells. Modules in the battery 
pack are subjected to equal amount of cooling air as shown in 
Fig. 2. Hence, it is sufficient to study the thermal response of a sin¬ 
gle module. SFTP driving cycle is more aggressive than Urban 
Dynamometer Driving Schedule (UDDS) and The Highway Fuel 
Economy Test (HWFET) and more energy is being charged to and 
discharged from the battery pack and more heat is generated. 
Therefore, it was chosen for analyzing the thermal response of 
the battery pack and benchmark with 18,650 cell battery pack. 
Thermostat feature was included in the battery model with switch¬ 
ing on the temperature set at 35 °C. After the cell temperature 
reaches 35 °C, the blower will come into deliver 10, 25 or 100 
CFM of cooling air per module. The specific parameters of the EV 
used for driving cycle simulations are tabulated in Table 1. 

2.4. Experimental setup and battery model parameters extraction 

Two different sizes of commercial cylindrical LFP cell were 
examined in the present study. The details of the cells are given 
in Table 2. Battery cycler (Maccor Instrument 4000) was used to 
charge and discharge the battery to obtain the voltage and current 
curve in function of time. 2.3 V was used as the cutoff voltage for 
constant current discharging while charging of the batteries were 
done in two modes: constant current charging until 3.65 V at 0.1 
l t followed by constant voltage charging until the current dropped 
to 0.01 A. Slow charging is desired to ensure the chemical process 
within the cell occurs at a rate similar rate for the transfer of elec¬ 
tric energy in the isothermal state. The batteries were put at rest 
for 1 h to ensure the temperature of the battery is in equilibrium 
with room temperature and the open circuit voltage of the cell is 
stable. Impedance analyzer (Solartron analytical 1400) was used 
to assess the internal resistance (R int ) of the batteries. The param¬ 
eters for 38,120 and 18,650 battery models are provided in Table 2. 

The surface temperature of the batteries was measured with 
twelve thermocouples (T-type) attached to different locations of 
the cell. Three thermocouples were attached in the axial direction 
and four sides of the battery surface. The experiments were con¬ 
ducted at room temperature of 25 °C under natural convection. 
Temperature readings were recorded using the HP34970A data 
acquisition system. The specific heat capacity of the cell was mea¬ 
sured using adiabatic accelerating rate calorimeter (THT ARC) as 
shown in Table 2. 

2.5. Numerical procedures 


Qgen — i[ Eq — V b att + T 


i 2 R c 


( 6 ) 


The electro-thermal model of the battery was developed using 
Matlab-Simulink 2011b. Electrical and thermal response of the 
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Exhaust Air 
Outlet: T air 


Cross Section 
for internal flow 


Fig. 2. Schematic of active air cooling system of an EV battery pack [17]. 


Table 1 

Vehicle and cooling system specific parameters for 38,120 cell. 


Parameter 

Value 

Vehicle mass, kg 

1828 

Frontal area, m 2 

3.238 

Coefficient of drag, Cd 

0.35 

Electric motor 

75 kW, 200 Nm max 

Battery pack 

19.5 kWh 

Number of cell per module 

28 

Number of module 

28 

A m f per module, m 2 

0.00868 

A ms per module, m 2 

0.401 

Mass flow rate of cooling air for battery pack, cfm 

280, 700, 2800 

Surrounding temperature, °C 

30 


Table 2 

Details of LFP cells used in this study. 

Parameters 

Cell 



18,650 

38,120 

Nominal voltage, V 

3.2 

3.2 

Nominal capacity, mAh 

1300 

8000 

Cathode material 

LiFeP0 4 

LiFeP0 4 

Anode material 

Graphite 

Graphite 

Terminal connector 

Spring loaded 

Screw 

Diameter, m 

0.018 

0.038 

Length, m 

0.065 

0.0146 

Weight, kg 

0.030 

0.355 

Specific heat, J kg -1 K 1 

900 

998 

E 0 , V 

3.21 

3.27 

R, Q 

0.03 

0.0034 

K, Q or V(Ah)- 1 

0.0119 

0.00216 

A, V 

0.2711 

0.0854 

B, (Ah)” 1 

152.130 

23.097 


battery model was extracted and compared with the experimental 
results of constant current discharging and Simplified Federal 
Urban Driving Schedule (SFUDS) dynamic power profile. Next the 
validated battery model is further extended to investigate the 
thermal response of the battery pack for a converted EV under 
US06-SFTP driving cycle. 

2.6. Cell selection for EV battery pack 

The electro-thermal model results of 38,120 cell and 18,650 cell 
will be used to determine the suitable cell for EV battery pack. The 
benchmarking process is based on Hyundai Trajet and thermal 
response of the battery pack is assessed the electro-thermal model 
under US06-SFTP driving cycle. The details of the results are 
provided in Section 3.6. 

3. Results and discussion 

3.1. Validation of the cell voltage 

Calibration result of the 38,120 battery model using 0.2 I t of 
discharge curve is shown in Fig. 3(a) while the error of prediction 
is given in Fig. 3(b). The highest error of the simulation results as 



Fig. 3. (a) 38,120 cell voltage during discharge at 0.2 I t to extract the model 
parameters, (b) Error of voltage prediction. 


□ II t Experiment 



compared to experimental data is about -0.06 V during initial 
discharge, while at the end of discharge it is less than 3%. In 
Fig. 4 the measured and simulated voltage for 38,120 cell during 
1, 2 and 3 I t -rate of discharge is shown. Comparisons of simulated 
and measured voltage show that the battery model provides a good 
estimation of the electrical behavior of the cell at various J r rates. 
Average squared error for 1, 2, 3 J r rates are 0.0000260, 
0.0000419 and 0.000117 respectively. Discharging at high J r rates 
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could cause a rise in cell temperature. In addition, internal 
resistance of the cell is dependent on temperature, thus the error 
of prediction at 3 I t is slightly larger as compared to 1 and 2 J r rates. 
Although the cell nominal rating is 8 Ah, the final capacity resulted 
from 1 I t of discharge could reach 8.6 Ah. The final capacity at the 
end of 3 J t of discharge is about 8.4 Ah which is about 3.5% less than 
the capacity at 0.2 I t (8.7 Ah). On the other hand, at 3 l t discharge of 
18,650 cell suffered more capacity loss which is about 8% as 
compared to its capacity of 0.2 I t [32]. 

3.2. Evolutions of the cell temperature and heat generation 

Simulated and measured average skin temperature of 38,120 
cell at different J r rates of discharging under natural convection 
(5 W m -2 K -1 ) is shown in Fig. 5. The results showed that at high 
I r rates of discharging, temperature rise of the cell are significant. 
The skin temperature of the cell could rise to 13.8 °C at 3 I r rate 
and the cell exhausted its energy in 20 min. Natural convection is 
no longer sufficient to keep the cell temperature within safe oper¬ 
ating limits. Therefore, a battery thermal management system is 
required in this situation. The averaged squared error of the 
experiments data and simulation results for 1, 2 and 3 J r rates of 
discharging are 0.0888, 0.0224 and 0.0143 respectively. At 3 I t of 
discharge, the skin temperature of the 38,120 cell is about 6 °C less 
than 18,650 cell at a similar rate of discharge. Heat exchange sur¬ 
face area per unit jelly roll volume of the cell is a critical factor to 
determine the heat dissipation rate from the cell to the environ¬ 
ment [39]. The heat transfer area per unit volume of 18,650 and 
38,120 cells is 222.22 m _1 and 105 m -1 respectively. As the ratio 
of surface area to volume decrease, the heat transfer ability of 
the cell was reduced and the internal temperature of the cell is 
increased [9]. Although the skin temperature of the 38,120 cell is 
lower than 18,650 cell, more heat is trapped inside the cell. There¬ 
fore, the smaller cell will have a higher skin temperature at a sim¬ 
ilar I r rate of charging and discharging. Sometimes, this will give 
false information to the user that a large cell performs better than 
a smaller cell in the thermal aspect. 

Estimated total heat generation of 38,120 battery is depicted in 
Fig. 6. 6.5 W, 14 W and 20 W per cell of heat were generated 
toward the end of 1-3 J r rate of discharge respectively. Heat gener¬ 
ated from the battery is also positively correlated to the size and 
capacity of the cell. Heat generated from 18,650 cell is only 9% of 
38,120 cell at 1 I t of discharge, while the heat generated at 3 I t is 
about 15% of 38,120 cell in a similar J r rate of discharge [28]. 

3.3. Internal temperature of cell 

Temperature distribution of 38,120 and 18,650 cells under 
different cooling conditions at the end of 3 I t of discharge is 




Fig. 6. Heat generated of 38,120 cell during discharge at 1, 2 and 3 I t - rates. 


provided in Fig. 7. Due to the anisotropic nature of layered active 
material inside the battery, the temperature distribution in the 
axial direction of the cell is fairly uniform [39] and a temperature 
gradient is expected only in the radial direction of the cell. Heat 
shrink film which is used for insulation of the metal casing will 
prevent the heat generated in the cell from being effectively dissi¬ 
pated to the environment. Therefore, a temperature jump was 
found between the insulator and metal casing of the cell as shown 
in the initial and end section as in Fig. 7. The maximum tempera¬ 
ture region of the cell is located in the active material region near 
the hollow core. The hollow core situated at the center of the cell is 
a product of the cell manufacturing process. During cell sealing 
process, an electrolyte will be injected into the active material 
and the hollow core will be filled with the electrolyte as well. An 
attempt to measure the internal temperature of the cell at the hol¬ 
low core region may not produce any useful information on the 
maximum temperature of the cell [38,39]. As shown in Fig. 7, the 
temperature in the hollow core region is lower than the active 
material near the hollow core region. 5 °C of temperature gradient 
was developed in 38,120 cell during 3 J t of constant current 
discharge under natural convection as shown in Fig. 7(a). The max¬ 
imum internal temperature of the cell may reach 44 °C for 3 I t of 
constant current discharge. On the other hand, the temperature 
gradient in the 18,650 cell is also 5 °C under natural convection 
(5 W m -2 K 1 ) and the maximum temperature in the active mate¬ 
rial region is about 50 °C. The temperature gradient of the cell will 
grow when strong forced convection is used to cool. When the heat 
transfer coefficient increase 500 W m -2 K 1 , the temperature gra¬ 
dient for 38,120 and 18,650 cell may rise to 14 °C and 16 °C respec¬ 
tively as shown in Fig. 7. Opposite to excellent thermal 
conductivity of metal casing, the active material region is a poor 
thermal conductor. Although the current collectors of the cell are 
made of copper and aluminum, but porous electrodes and separa¬ 
tors are poor conductors. The main function of a separator is to iso¬ 
late the anode and the cathode layer while providing an effective 
transport medium for Li-ion. Poor thermal properties of the sepa¬ 
rator will prevent the heat from being effectively dissipated to 
the outer environment and it further proves that heat flux inside 
the large cell is not always outward [39]. Therefore, the safety of 
the battery cannot be simply determined by measuring the skin 
temperature of the cell. Using strong forced convection to cool 
the cell is not encouraged as this will introduce undesirable tem¬ 
perature gradient in the cell and accelerate the thermal aging. In 
the open literature, there are various types of degradation models 
used to investigate the cycle life of the battery [40-42]. In the cur¬ 
rent study, thermal aging model proposed by Wang et al. is 
adopted and the thermal aging of the LiFeP0 4 cell can be repre¬ 
sented by the Arrhenius equation as below [40] : 

(loss = 30330 exp (^|^)< 552 (9) 
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Fig. 7. Predicted variation of internal temperature of the cell, (a) 38,120 cell, (b) 18,650 cell at 3 l r rates. 


The capacity loss of the cell at 30 °C, 40 °C, 50 °C and 60 °C for 
total of 1000 Ah throughput is 5.12%, 7.65%, 11.12% and 15.81% 
respectively. This further indicates that 10°C of temperature 
difference will lead to additional 3-5% reduction in the cycle life 
of the cell. Compared to strong convection cooling of the battery 
outer surface, improving the thermal properties of the separator 
and internal structure of the cell or filling the hollow core with 
heat absorbing medium is more effective in improving the heat 
dissipation of the cell. 

3.4. Dynamic behavior of the cell under SFUDS 

Validation of battery model dynamic behavior required com¬ 
parison with an independent set of series charging and discharging 
test. SFUDS, featuring 360s of repeating charging and discharging 
at certain specific power was used for this study [28]. The voltage 
and current response plot of the simulated results and the 



experimental measurement of 38,120 cell in the first 360s are 
depicted in Fig. 8. The error graph shows that the discrepancy 
between the model and experiment for the voltage is below 0.5% 
and below 2% for the current. Fig. 9 shows both experimental 
and simulated voltage and temperature profile of 38,120 cell under 
450 min of SFUDS cycle. The voltage drops to the lowest point dur¬ 
ing 79W/kg specific power discharge within each sub-cycle of 
360s. The voltage of the cell decreases with the number of sub¬ 
cycles and oscillates in a pulse with the power pulses. The voltage 
of the cell gradually decreases as the SOC of the cell ceases to cut 
off voltage to maintain the power requirement of the cycles. A good 
agreement can be observed between the simulated and experi¬ 
mental results. Total heat generated for 450 min of a SFUDS test 
of 38,120 cell is about 106J which can be effectively dissipated 
by natural convection and the overall change of the cell skin tem¬ 
perature is within 1 °C. Endothermic heat generation is observed 
during the initial 280 min of the cycle. Then, the heat generation 



Fig. 8. Validation of battery model using SFUDS (a) Voltage, and (b) current. 
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Fig. 9. Comparison of experiment and modeling results for 450 min of SFUDS 
profile, (a) Voltage, and (b) temperature. 


rate gradually increased and became exothermic as the charging 
and discharging process continued, contributing to the rise of the 
overall cell temperature. Therefore, the battery model is crucial 
to provide an insight of battery state during testing in an economic 
and nondestructive way without over-charging or over¬ 
discharging of the cell. 

3.5. Thermal response of the battery pack under US06 driving cycle 

Thermal response of the 38,120 cell battery pack under US06 
driving cycle is shown in Fig. 10. The battery pack is subjected to 
600s of the repeated US06 cycle until the voltage of the battery 
pack reached 64 V. 38,120 cell battery pack is able to complete 
4.5 cycles of US06 and deliver power for 45.1 min. Temperature 


of cells increase continuously in the cycling process and reaches 
the maximum at the end of the cycle. At the end of the cycle, the 
average temperature of the battery in the module is 38.1 °C, 
37.6 °C and 36 °C for 10 CFM, 25 CFM and 100 CFM of air flow 
per module respectively. The total flow rate of air is 280, 700 or 
2800 CFM for 38,120 cell battery pack which consists of 28 mod¬ 
ules. As compared to 18,650 cell battery pack in the previous study 
[28], 38,120 cell battery pack has higher energy used, higher 
energy gained from regenerative braking and higher capacity. 
The heat generated from 38,120 cell module is not as intensive 
as compared to 18,650 cell module in the previous study [28]. As 
shown in Table 3, heat generated from 38,120 cell battery pack is 
about 33% less than the 18,650 cell battery pack. This is caused 
by additional heat generated from the contact resistance of spring 
loaded terminal connector for 18,650 cells. On the contrary, screw 
connector is used for 38,120 cells to connect the cell and the con¬ 
tact resistance is very low (~0.6 miliohm), only 0.45 kj of heat is 
generated by the contact resistance throughout the cycle for 
38,120 cell battery pack. On the other hand, the effect of contact 
resistance of 18,650 cell pack contributed 0.764 MJ of heat 
generated. 

At initial 2200s, natural convection is sufficient to dissipate the 
heat generated from the 38,120 cell battery pack in the US06 cycle. 
The blower came in after 2200s (3.6 cycles) to supply cooled air to 
chill the cells and attempt to bring down the cell skin temperature 
to below 35 °C. Conversely, skin temperature of 18,650 cell reached 
35 °C as soon as 1.2 cycles (744s) of the US06 cycle were 
completed. Cross section flow area for cooling air per 38,120 cell 
module is 0.00868 m 2 and is about twice the flow area of 18,650 
cell module. Hence, mass flow rate of cooling air for 38,120 cell 
module need to be double to achieve a similar cooling condition 
as 18,650 cell module. Although the blower is only operated for 
8.3 min, large mass flow of air (2800 CFM per pack) is needed to 
suppress the temperature rise of the cell. In this case, huge blower 
is needed and this will increase the cost and power consumption of 
the pack thus reducing the driving range of the EV. 

Gradient of internal cell temperature under strong forced con¬ 
vection, 25 CFM for 18,650 cell module and 100 CFM for 38,120 cell 
module was investigated. As shown in Fig. 11, gradient of the 
38,120 cell internal temperature could reach 6°C, although there 


Table 3 

Energy distribution for US06 cycle per pack. 


Cell 

Discharging, MJ 

Charging, MJ 

Heat generation, MJ 

18,650 

58.99 

5.71 

4.937 

38,120 

75.73 

7.40 

3.319 



Fig. 10. Thermal response of 38,120 battery pack under US06 driving cycle. 
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Fig. 11. Gradient of cell internal temperature for US06 driving cycle under 
maximum cooling capacity. 



Fig. 12. Cell skin temperature and maximum internal temperature of 38,120 cell 
under laminar flow of ethylene glycol. 

is only 4 °C of temperature difference was developed in the 18,650 
cells at the end of the US06 cycle. The internal temperature of 
18,650 cell and 38,120 cell could reach 38 °C and 39 °C respec¬ 
tively. Hence, a cell with a larger diameter will have a larger 
temperature gradient as compared to a small cell. 

The effect of using liquid cooling on the 38,120 cell battery pack 
was also investigated. For closed pack architecture of 38,120 cell 
module, laminar flow of ethylene glycol is sufficient to maintain 
the temperature rise of the cell surface below 5 °C as illustrated 
in Fig. 12. Although the skin temperature of 38,120 cell is success¬ 
fully maintained at 33 °C, large temperature gradients develop 
inside the cell. The internal temperature of the cell could reach 
42 °C which brings adverse effect to the cell cycle life and acceler¬ 
ates the thermal aging of the cell. 

On the other hand, immersion cooling using 3 M Novec 7000 
engineering fluid or Fluorinert electronic liquid could be a solution 
for thermal management of large Li-ion cell. The cells are 
immersed in the liquid and take advantage of large latent heat of 
evaporation of the liquid to achieve cooling purpose. 

3.6. Cell selection and integration issues ofEV battery pack 

The selection and integration of cells to form a battery pack are 
important to determine the safety, durability, life cycle, driving 
range, cost and fuel economy of the EVs or HEVs. Poor integration 
of cells may lead to disadvantages such as low performance of the 
battery pack, shorter cycle life, safety hazard and unnecessary pro¬ 
duction costs [43,44]. There are several issues associated with the 
integration of cells into a battery pack, and this involves the 


mechanical, electrical, thermal, control and monitoring, packaging, 
manufacturing, maintenance and repair aspects. These issues are 
important to create a comfortable environment to bring out the 
best of the cells. Some automakers prefer using a large number 
of small cells for battery pack while others prefer using a few large 
cells. For example, Tesla Roadster used more than 7000 18,650 
cells (3100 mAh) for their Tesla Model S while Mitsubishi used 
88 much larger cells (50 Ah) for their Mitsubishi i-MiEV. Using a 
larger number of small cells or a small number of large cells has 
its own advantages and disadvantages. The advantages of using 
small cells include cost efficiency, a lack of thermal aging tendency 
(more uniform temperature) and improvements in safety. Con¬ 
versely, the disadvantages include many interconnections, higher 
integration and assembly cost, lower weight and volume efficiency 
and lower reliability. In addition, a large quantity of cell protection 
and monitoring circuit boards are needed. [43]. 

On the other hand, using larger cells has several advantages 
such as lower assembly cost, higher weight and volume efficiency, 
higher reliability, less complex interconnections and ease of trou¬ 
bleshooting [43]. The disadvantages include higher cell production 
costs, a tendency to thermal aging (high gradient of temperature 
across the cell), lower quality and increased difficulties in thermal 
management. The system will also be prone to safety issues and 
also needs a complicated and expensive battery protection circuit 
(high ampere current sensor and used of MOSFET instead of relays) 
[43]. 

In this section, integration issues of two different sizes of 
cylindrical cell will also be explored. The benchmarking of the 
38,120 and 18,650 cells is based on the converted EV using the 
Hyundai Trajet. A relative score of “better than” (+) or “worsen 
than” (-) is used to represent how the subjective criterion is rated 
in comparing two different sizes of cells. The assembly of 38,120 
and 18,650 cell unit is shown in Fig. 13. The comparison of the cell 
integration issues from various aspects is tabulated in Table 4. 

Manufacturing technology for small cells is more mature than 
the large cell and the price is cheaper for mass purchasing. In addi¬ 
tion, the total number of 18,650 cells (5040 cells) is 6 times higher 
than the total number of 38,120 cells (784 cells) in the designed 
battery pack. The total weight of 18,650 cells is 201.6 kg, while 
the total weight of 38,120 cells is 278.32 kg. This is due to the thick 
metal pieces used for 38,120 cell threaded terminals. The total 


Fig. 13. Assembly of the cells. Left: 18,650 cell. Right: 38,120 cell. 
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Table 4 

Selection and integration issues of 18,650 and 38,120 cell for battery pack. 


Criteria 

Cell 

18,650 

38,120 

Packaging 

Number of cell 

5040 

784 

Weight, kg 

201.6 

278.32 

Volume, m 3 (closed pack) 

0.106 

0.150 

Weight of bus bar, kg 

20.26 

12.45 

Weight of cell holder, kg 

-30 

-9 

Weight of terminal connector, kg 

12.13 

- 

Cell cost, USD 

3-5 

-20 

Assembly of single cell 

a, (3 orientation of cell [36] 

360°, 0° 

360°, 0° 

Cell handling + insertion time, s [36] 

3.5 

3.5 

a, (3 orientation of bus bar [36] 

180°, 

180°, 


180° 

180° 

Bus bar handling + insertion time per cell, s [36] 

15.72 

15.72 

a, [3 orientation of cell holder [36] 

360°, 

360°, 


360° 

360° 

Cell holder handling + insertion time, s [36] 

7.4 

7.4 

a, [3 orientation of terminal connector [36] 

360°, 0° 

- 

Terminal connectors handling + insertion time, s [36] 

6.6 

- 

Bus bar assembly time (two terminal), s [36] 

46.36 

46.36 

Assembly cost, USD (assumed 5 USD per hour) 

0.111 

0.101 

Electrical and control 

Terminal contact resistance, mQ 

10 

0.6 

Complexity of wiring 

1 

7/45 

Cell monitoring 

- 

+ 

Reliability 

+ 

- 

Thermal management 

Gradient of internal temperature, °C (Based on US06 
cycle) 

4 

6 

Heat generated from contact resistance, kj/cycle (based 
on US06 cycle) 

208.674 

74.08 

Heat generated from the battery pack, kj/cycle (based 
on US06 cycle) 

870.485 

780.543 

Power consumption for cooling fan 

1 

0.77 

Complexity of cooling system design, cell/m 2 

3086 

693 

Maintenance and troubleshooting 

Time required to identify the faulty cell in the parallel 
bank 

1 

7/45 

Ease of cell replacement and services 

1 

7/45 


weight of 18,650 cell battery pack increased to 264 kg when bus 
bar, terminal connectors and plastic cell holder are installed, while 
features such as the battery pack support structure, thermal 
management system, battery protection circuit board and wiring 
are excluded. The total weight of the 38,120 cell battery pack is 
about 300 kg and it is much heavier than 18,650 cell battery pack. 
Furthermore, smaller cells have better packing density and power 
density as compared to larger cells. Therefore, from the packaging 
point of view, smaller cell is better for battery pack integration. The 
weight of large cell could be improved by using different material 
for battery casing and cell terminals. 

Most of the products are manually assembled and Boothroyd 
Dewhurst DFA method for manual assembly is widely used to 
assess the assembly efficiency. More than 50% of the total manu¬ 
facturing cost and nearly 60% of the total production time are spent 
on the assembly process [45-47]. For manual assembly study, the 
analysis is split into two distinct analyses, which are known as 
manual handling and insertion analysis respectively. In this study, 
orientation, part size, ease of handling, types of cell fastener, bus 
bar and terminal connector are examined to determine the 
handling and insertion time, a rotation refers to the alignment of 
the axis of the part corresponded to the axis of insertion, while 
the rotation of the part about its axis of rotation is called p rotation. 
In the analysis, the assembly of a single cell is taken into consider¬ 
ation. The labor cost of the assembly is assumed as USD 5 per hour. 


The total assembly cost of the battery pack excluding battery man¬ 
agement system can be calculated by multiplying the number of 
cell, battery holder and bus bar used. The analysis results are 
tabulated into the assembly section of Table 4. The total assembly 
cost of a single 18,650 cell required USD 0.111 while the assembly 
of a single 38,120 cell required USD 0.101. For the 18,650 cell 
battery pack that consists of 5040 cells, about USD 559 and 
112 h of workmanship are needed for cell assembly and this 
excludes the assembly of thermal management system, battery 
protection circuit board and wiring. On the other hand, a 38,120 
cell battery pack only required USD 79 and 15.84 h of workman¬ 
ship. Despite the differences in assembly time, the 18,650 and 
38,120 cells actually share the similar characteristic of assembly 
in term of a and p rotation. 

Another difference in the mechanical aspect is the use of 
terminal connectors. Threaded terminals are provided for 38,120 
cells and no terminal connectors are required to assemble 38,120 
cells. On the other hand, the 18,650 cells require spring- loaded 
terminal connectors to bridge the cell terminal to the bus bar. 
Hence, an additional component and assembly process are 
involved for the 18,650 cells. This will incur extra manufacturing 
and assembly cost. As such, in the assembly point of view, the large 
cell has better assembly efficiency than the smaller cell and the 
total production time can be shortened effectively to produce a 
more economical battery pack. 

From the electrical aspect, an increasing number of cells also 
mean that the number of terminal connectors, bus bars and the 
interconnections are proportionally increased. Furthermore, this 
will increase the complexity in wiring, copper loss, causing 
difficulty in cell monitoring and protection circuit design. In the 
consideration of design for assembly and disassembly, spring- 
loaded terminal connectors are chosen over spot welding on nickel 
strip for 18,650 cells assembly. Nevertheless, it will produce a 
higher value of contact resistance and more energy is wasted in 
copper loss, which means that the total usable energy of the bat¬ 
tery pack is reduced. Moreover, a large number of cells also indi¬ 
cate a large quantity of cell protection and monitoring circuit 
board is needed. Although a large cell requires less cell protection 
and monitoring of the circuit board, the circuit design is more com¬ 
plicated and involves high current electronic devices such as the 
high ampere current sensor and thick interconnections. Moreover, 
cell balancing (active and passive) is more difficult for a large cell. 
From the electrical point of view, wiring for a large cell is less com¬ 
plex and copper loss can be reduced to a minimum. Furthermore, 
the large number of interconnections will increase the chances of 
failure and make it more difficult for troubleshooting. Hence, from 
the control point of view, fewer cell protections and monitoring 
board are needed for large cells and this makes the control task 
easier. On the other hand, having more cells connected in parallel 
is expected to increase the reliability of the battery pack and the 
effect of bad cell could be minimized [48,49]. Cells with the same 
chemistry have the same mean time between failures. For exam¬ 
ple, either size of the cells has about 1% likelihood of sudden lose 
of capacity due to internal failure or external wiring problem over 
a month. The studied showed that after 10 months, the battery 
pack with more cells connected in parallel is less influenced by 
capacity fading [48,49]. 

US06 driving cycle is used as a reference to assess the thermal 
performance of 38,120 and 18,650 cell battery packs. The heat gen¬ 
erated from 18,650 and 38,120 cells is 870.485 kj/cycle and 
780.543 kj/cycle respectively. In the 18,650 cell battery pack, the 
heat is shared equally in 5040 cells. On the other hand, in the 
38,120 cell battery pack, the heat is only shared with 784 cells 
and the thermal problem is more serious for large cell. Due to 
the small number of interconnections in the 38,120 battery pack, 
the heat generated from the contact resistance is also small which 






106 


L.H. Saw et al./Applied Energy 131 (2014) 97-107 


is about 74.08 kj/cycle and 2.8 times less than the heat generated 
from the contact resistance of the 18,650 cell battery pack. Power 
consumption for cooling fan is calculated using Eq. (10) [ 17 and an 
18,650 cell battery pack is used as a reference for comparison. 

Vced = 30(^) 08 ( 10 ) 

Axial cooling architecture is used to benchmark the power con¬ 
sumption for the cooling fan. 18,650 cells are stacked in two layers 
to achieve the similar length with the 38,120 cells. Due to poor 
packing density of the 38,120 cell battery pack as discussed in 
Section 3.5, for a given mass flow rate of cooling air, heat transfer 
coefficient is only 0.77 times less than the 18,650 cell battery pack. 
Conversely, using a large number of small cells rapidly increases 
the complexity and the cost of the battery thermal management 
system design and the assembly of the cells [43]. For the axial cool¬ 
ing architecture, there will be 3086 pieces of 18,650 cell and 693 
pieces of 38,120 cells in one square meter. The manifold of the 
18,650 cell battery pack needs to take care about 4.5 times more 
cells than the manifold of the 38,120 cell battery pack. Hence, 
the design of the 18,650 cell battery pack manifold and ducting 
is more complex and challenging to ensure the uniformity of tem¬ 
perature distribution. The only concern is that the thermal man¬ 
agement for a large cell is the thermal resistance in the cell. Heat 
is trapped in the cell and this leads to the development of a large 
temperature gradient. This produces an uncontrollable cell internal 
temperature, which will give rise to safety and thermal aging 
issues. Hence, from the thermal perspective, a cell with a large sur¬ 
face area to volume of active material is desired, such as thick and 
short cylindrical cell or thin and long cylindrical cell as it could 
help to dissipate the heat more effectively. On the other hand, long 
electron current path of the thin cylindrical cell, may lead to per¬ 
formance degradation and shorten the cycle life of the cell [43]. 
The placement of the current tab in thick and short cylindrical cells 
is a critical issue and needs to be optimized to result in an even 
current distribution in the active material. Besides, making the cell 
thinner, such as by changing to the pouch cell design, is more 
favorable for thermal management. 

In addition, a large number of cells make the troubleshooting 
and maintenance jobs of the battery pack tedious. This type of bat¬ 
tery pack is more suitable for use as in “use and throw” products. It 
does not encourage recycling or repair, and is not environment 
friendly. On the other hand, the replacement of faulty cells for 
the battery pack using large cells and the servicing of faulty cell 
is also easier. For example, in the 18,650 cell battery pack, there 
will be 180 cells connected in parallel to form a module. On the 
other hand, there will be only 28 cells connected in parallel for 
the 38,120 cell battery pack. The time required to detect the faulty 
cell for the 18,650 cell module is 6 times more than 38,120 cell 
module. Hence, a large cell is more favorable for maintenance point 
of view. 

By reviewing the integration issues of 18,650 and 38,120 cells, it 
is apparent that the large cell has extra advantages over the smal¬ 
ler cell. It should also be noted that the large cells tend to develop, 
a large gradient of temperature, which is associated with thermal 
aging issues. This problem, however, could be easily alleviated 
using an internal cooling approach. 

4. Conclusions 

A detailed battery model is developed to investigate the perfor¬ 
mance and thermal response of two different sizes of the cylindri¬ 
cal cell. The model is based on the modified Shepherd model by 
extracting the data obtained from 0.2 / r rate of constant current 
discharge curve is able to predict dynamic behavior of the cell with 


good accuracy. Detailed information about the battery operating 
parameters such as SOC, I-V characteristics, skin temperature 
and internal temperature of the cell can be obtained. The simula¬ 
tion results showed a good agreement with experimental data 
under various operating conditions. Heat generated in the cell is 
positively correlated with the J r rates and diameter of the cell. Heat 
exchange surface area per unit volume of the cell is a critical factor 
that determined the heat dissipation rate from the cell to the exter¬ 
nal environment. 

As the ratio of surface area to volume is reduced, the heat trans¬ 
fer decreases and the internal temperature rises. The maximum 
temperature region inside the cell is located in the circular region 
of active material near the hollow core. Due to the large thermal 
resistance and insulation effect of the separator, the temperature 
difference of the cell in the radial direction is significant and it 
increases with the diameter of the cell. This contributes to the slow 
rise of large cell skin temperature. Smaller sizes have better 
temperature distribution in the cell. Therefore, measuring skin 
temperature of the cell is not sufficient as a safety reference for cell 
operation. 

Strong forced convection cooling should be avoided as it will 
increase the cell internal temperature gradient and accelerate the 
rate of thermal aging. Other than cost, the electrical, thermal, 
mechanical, cell assembly and maintenance aspects are the issues 
that need to be addressed in the cells integration. These issues will 
further determine the performance, cost and life cycle of the bat¬ 
tery pack. Smaller cells allow easier thermal management, but 
result in a high number of interconnections and an increase in 
the chance of failure. Large cells perform better in term of assem¬ 
bly efficiency, cell monitoring and maintenance, but perform 
poorly in thermal management. Due to the growing in Li-ion cell 
industry and manufacturing technology, cell cost is no longer an 
issue. An improvement of the cell internal structure, however, is 
still a need to allow the heat to be dissipated quickly so that a large 
size and capacity cell can be built. 
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